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INTRODUCTION 


“Roughly, by a complex system I mean one made up of a large number of parts that interact 
in a nonsimple way” is Simon’s description (ref. 1) of a complex system. He goes on further to say 
that “complexity frequently takes the form of hierarchy,” and by intuitive arguments he points out 
that the evolution of complex systems is highly reliable if it is carried out as a hierarchic process 
whereby complex systems are formed by interconnecting stable simple parts (subsystems). 

The main objective of this work is to show rigorously that a complex system with (or designed 
to have) a competitive structure has highly reliable stability properties. The competitive models were 
studied in such diverse fields as economics (refs. 2, 3) and biology (ref. 4), arms race (ref. 5), and 
pharmacokinetics (ref. 6), and only recently it was shown (refs. 7—10) that in these various scientific 
disciplines, such models are “fail-safe” stable. In the framework of the connective stability concept 
(refs. 11-14), we will provide a definite support of Simon’s intuitive arguments. We will show that, 
under certain conditions, a stable complex system can tolerate a wide range of nonlinear, time- 
varying perturbations. 

For large-scale systems, a competitive hierarchic model will be constructed by aggregating the 
stability properties of each subsystem so that stability of the model implies stability of the original 
system despite structural perturbations whereby subsystems are disconnected and again connected 
in various ways during the operation of the system. That is, from stability of each subsystem and 
stability of the aggregate model on the upper hierarchic level, we infer connective stability of the 
overall complex system. This result is remarkable in that it provides a natural setting for designing 
large-scale dynamic systems as competitive structures with highly reliable stability properties. 

The stability investigations of competitive models are carried out by using the powerful mathe- 
matical concept of the comparison principle and vector Liapunov functions (ref. 15). The concept 
is extended here to include considerations of the connective stability aspects and the effects of the 
structural perturbations. Furthermore, the decomposition-aggregation methods developed in this 
context can take advantage of the special structural features of complex systems and reduce con- 
siderably the dimensionality of relevant stability problems. 


♦The research reported here was supported by NASA Ames Research Center, Moffett Field, Calif., under 
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COMPETITIVE MODELS 


Let us start with a linear constant dynamic system described by the differential equation* 

x=Ax (1) 

where x(t) G R n is the state of the system and A = (dy) is an n X n constant matrix. By use of the 
connective stability concept (ref. 12), we derive the conditions under which stability is a highly 
reliable property of the system ( 1 ), and show that it remains stable despite a wide range of nonlinear, 
time-varying perturbations. In fact, by use of the modem mathematical machinery of the compari- 
son principle (ref. 15), we show that, under relatively simple conditions, stability of the system (1) 
implies stability of a broad class of dynamic systems described by 

x = A(t, x)x . (2) 

In equation (2), x(t) E.R n is again the state of the system and the n X n matrix function A: 

T X R n -*R n is defined, bounded, and continuous on T X R n so that the solutions x(/; t 0 ,x 0 ) of 
equation (2) exist for all initial conditions (t 0 , x 0 ) G T X R n and t G T 0 . The symbol T represents 
the time interval (r, + °°), where r is a number or the symbol -°°, and T 0 is the semiinfinite time 
interval (/ 0 , + °°). The matrix A = (a^) of equation (2) is an obvious modification of the matrix 
A = (dj.) of equation (1), in which the constant elements dy of A are replaced by nonlinear, time- 
dependent functions ay = ay(t, x). 

To consider the connective aspect of stability, the elements ay of the matrix A in equation (1) 
are written: 


x) = Syipft, x) + ey(t) \j/y(t, x) (3) 

where 8y is the Kronecker symbol and i pft, x), tyy{t, x) G C^° >°\t X R n ). In equation (3), 
e ij ~ e if t ) are e l ements of the n X n interconnection matrix E = (e^), which are ey(t) G C°(T) and 
are res/ricted as ey(t) G [0, 1 ] , Vt G T (ref. 12). For system (2), the element (t) reflects the 
coupling between x^t) and Xj(t) at each instant in time, that is, the time-dependent influence of the 
state xj(t) on the derivative of the state x^t). Therefore, the interconnection matrix E represents 
the structural perturbations of the nonlinear matrix system (2). 

In this section, we study asymptotic stability properties of system (2) under structural pertur- 
bations. More precisely, we investigate stability formulated as: 

Definition 1. The equilibrium state x = 0 of system (2) is connectively asymp- 
totically stable in the large if and only if it is asymptotically stable in the large 
for all interconnection matrices E(t). 


*With some obvious exceptions, lower case roman letters denote vectors, capital Roman letters denote 
matrices, and Greek letters denote scalars. 
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Before we turn to the derivation of the conditions for the kind of stability expressed by 
definition 1, we need the notion of the fundamental interconnection matrix E (ref. 12). The 
matrix E is a time-invariant, interconnection matrix in which the elements take on binary values 
— 1 if the / th state Xj influences the ith time derivative x ^ of the state x and 0 if Xj has no influence 
on x j-. The matrix E is a binary matrix (ref. 16) that reflects the basic structure of the system. 
Therefore, any interconnection matrix E(t) is generated from E by replacing the unit elements of 
E by corresponding elements e^(/) of E(t). 

The conditions for connective stability are expressed in terms of E , but are valid for all E as 
required by definition 1 . This is an important qualitative result since we show stability of a class of 
nonlinear, time-varying systems by proving stability of one member of that class which is a time- 
constant linear system. 

To establish conditions for asymptotic connective stability, we assume that the elements 

*) of the matrix A(t, x) are specified by equation (3) where \pj(t, x ) and x) are bounded 
functions on T X R n and that there exist numbers oq > 0, ot^ > 0 such that 

'PjU, *)i*/i > <*/$/( I*/ 0, x)Xj < (XjftjdXjO, ViJ = 1,2 V(t, x)GTXR n (4) 

and OLj > Oijj- . In equation (4), + are comparison functions that belong to the class 

H : 0,(p) e CO (/?|), 0 f (O) = 0 and 0 f (p) < 0 f (p 2 ), V Pi . p 2 : 0 < p, <p 2 < + °° (ref. 17, 18). 

By A = (a~), we denote the n X n constant matrix with the coefficients 

“ij = ~^ij°Hj + ' e i) a ij ^ 

where the elements e- take the values 1 or 0 according to the matrix E. 

We prove the following: 

Theorem 1. The equilibrium state x — 0 of the system (2) is connectively 
asymptotically stable in the large if the n X n constant matrix A = (a^j) 
defined by equations ( 4) and ( 5) satisfies the conditions: 



«l 1 

aj 2 

«1 k 




. 

if 

a 2 1 

#2 2 

a 2 k 






a ki 

^k2 

a kk 


>0, Vk = 1,2, ... ,n. 


( 6 ) 


Proof : Let us consider the function v: R n -*■ R+ , 


v(x)= £ bj\Xj\ 
i= 1 


(7) 
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as a candidate for Liapunov’s function (ref. 19) for system (2) where b- t > 0, i = 1, 2, . . . , n are yet 
unspecified numbers. 


For v(x), we have the inequalities 

0/( ifcc II) < u(x) < (Pjji llxll) , Vit, r)6TX R n 
where 0/ and 0// e // are given as 

0 / 11 * 11 ) = b m \\x\\ , 07/01x11) = n x n b M ||x|| 
and b m = mjn bj and bj^ — max b- . 


( 8 ) 


(9) 


Since the derivative of |x ; (t)l need not exist at a point where Xj(t) — 0, it is necessary to calcu- 
late the right-hand derivative D + Uq(f)t with respect to equation (2) as proposed in reference 19. 

For this purpose, the functional o 2 is defined as 

I I , if Xj > 0, or if x 2 - = 0 and x 2 - > 0 \ 

0, if x 2 - = 0 and Xj = 0 ? (10) 

-1 , if x ; - < 0, or if x 2 - = 0 and x < 0 , / 

where x 2 - = x^t) £C l (T ). Then 

D + \x i (t)\=a i x i (t). (11) 

Using the constraints (4) and expression (1 1), we calculate the desired derivative as 


D + v(x) = ^2 bf-x^t) 
i= 1 


=1L b i°i X) 

i= 1 /=! 


x)xj(t) 


<b T Aw(x), V(t, x) & T X R n 


where b = (6, b 2 ... 6^)^ is a positive constant vector (6 > 0) and the positive vector function 


( 12 ) 


w: R n R + n is defined as 


w(x) = [ 0 ! ( |X ! |)0 2 (|x 2 |) . . . <p n (\x n \)] 


(13) 


From equation (5), we conclude that A is a Metzler matrix (ref. 20), that is, it has negative 
diagonal elements (a^ < 0) and nonnegative off-diagonal elements (a^ /). We recall that a 

matrix A is called a Hicks matrix (ref. 3) (equivalently, —A is an M matrix (ref. 21) if all even-order 
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principle minors of A are positive and all odd-order principle minors of A are negative. For a 
Metzler matrix A, the Hicksian property is equivalent to the Sevastyanov-Kotelyanskii conditions (6) 
(ref. 22). Since A is a Metzler matrix, the fact that A satisfies inequalities (6) and is a Hicks matrix 
is equivalent (refs. 20, 2 1) to saying that, for any constant vector c > 0, there exists a constant 
vector b > 0 such that 


J = -b T A. (14) 

Therefore, we can rewrite inequality ( 1 2) as 

D + i>(x) < -c^w(x) 
n 

05) 

1= 1 

<-0///(IWi), V(t,x)ETXR n 
where c m = min c ; - and 0///OMI) G//. 

From (8) and (15) and reference 18, we conclude global asymptotic stability of x = 0 in 
equation (2). To show that stability is also connective, we need only notice that 

A{t,x)x<Aw{x), V(t,x)ETXR n (16) 

where inequality (16) is taken component-wise. Therefore, equation (15) holds for all E(t). This 
proves theorem 1 . 

Note that the constraints (4) imply that i x ) > 0, V(t, x) ET X R n . Positivity of x) 
is absolutely essential for stability of equation (1) since it is easy to show that the Hicks condi- 
tions (6) imply < 0, Vi — 1,2 With this in mind, we can rewrite the first condition in 
(4) as |^(/, > <Xj<t>ji(\x^\), which looks similar to the second condition (4) except for the 

reversal of the inequality sign. 

If the conditions (4) are simplified to 

\p i (t,x)>a i , |i|/ I y(t,x)x / i yi,j = 1, 2, . . . , n\ V{t, x) E T X R n (17) 

where comparison functions <pj(\Xj\) are chosen as \Xj |, then we can establish exponential property 
of connective stability as in reference 9. Furthermore, if we use the notion of absolute stability for 
the nonlinear matrix systems proposed by Persidskii (ref. 23), we can prove that Sevastyanov- 
Kotelyanskii inequalities (6) become both necessary and sufficient conditions for connective 
stability (ref. 9). 

On the basis of the constraints (17), we define the following classes of continuous functions: 
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(18) 


i>i = : \p i (t,x)>a i \ 

► 

4>ij = : n//y(r,x)| <ofy-| 

- 

where Oj, o^- are numbers as in (4). Then, we state: 

Definition 2. The equilibrium state x = 0 of the system equation (2) is con- 
nective ly, absolutely , and exponentially stable if and only if there exist two 
positive numbers n and ir independent of initial conditions (t 0 , x 0 ) such that 

ll*(t;t 0 ,-*o)ll <n||x 0 ll exp[-7r(f-f 0 )] ,VtST 0 (19) 

for all (t 0 , x 0 ) E T X R n , all \J/j E and all interconnection 

matrices E(t). 


To establish this kind of stability, we can use the following: 

Theorem 2. The equilibrium state x — 0 of the system (2) is connectively, 
absolutely, and exponentially stable if and only if the n X n constant matrix 
A = (djj) defined by ( 5) and (17) satisfies conditions ( 6). 

Proof. Let us consider again the function !>(*) in equation (7). When (4) is reduced to ( 1 7), 
the vector w(x) in equation ( 1 3) becomes [ \x j | \x 3 \ . . . \x n ) ) T , and from ( 1 2) and ( 1 6) we get 

n n n 

D + v(x) < X bjlxjWdjjl + X Ixyl X , V(t, x)ETX R n . (20) 

7-1 /= 1 1=1 

&1 


Since A is a Metzler matrix, the fact that it satisfies conditions (6) is equivalent (refs. 21, 24) 
to saying that there exists a positive vector b = (b l b 3 . . . b n ) , and a positive number 7r such that 

n 

\djj\ ~ bp X b i$ij\ >«, Aj=l,2,...,n (21) 

z-1 

i*j 

that is, A is a quasidominant diagonal matrix (ref. 24). 


From (20) and (21), we get the differential inequality 

D + u<-ttv, VtET 0 , VvER + l , VE . (22) 

By integrating equation (22), we obtain 

^[*(01 <K*o) exp[— 7T(r — r 0 )] , VtET 0 , V(t 0 ,x 0 )E T X R n , VE . (23) 
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Using the well-known relationship between the Euclidean and absolute-value norms 
I hr 1 1 < |*| n ||;c|| , we can rewrite (23) as 


with 


||x(t; t 0 ,*o)ll <n||x 0 || exp[-7r(/-t 0 )] , Vt €E T 0 
V(t 0 ,x 0 )ETXR n , F^e* f , F^y G ^y , VE 

n=n i/2 b M b m - 1 


( 24 ) 


(25) 


where = max bj and b m = min b . 

Therefore, conditions (6) are sufficient for the absolute exponential property of connective 
stability of the equilibrium x = 0 in equation (2). This establishes the “if’ part of theorem 2. 

To prove the “only if’ part of theorem 2, we select the particular system (2) specified by 

x) = ay , \J/ t j{t , x ) = a t j , Vi, j = 1, 2, . . . , n (26) 

and the fundamental interconnection matrix E. That is, the matrix A(t, x) in equation (2) is taken 
as the constant Metzler matrix A and system (2) is described by equation (1). If A does not satisfy 
the Hicks conditions (6), the system (1) is unstable, and the equilibrium x = 0 of equation (2) is not 
stable F^E'Ej, F^y G 'Ey- . This completes the proof of theorem 2. 


HIERARCHIC MODELS 


On the basis of the results obtained in the preceding section, we conclude that stability is a 
highly reliable property of competitive dynamic systems since it can tolerate a wide range of struc- 
tural, nonlinear, and time-varying perturbations. Therefore, it would be desirable to define a class 
of noncompetitive dynamic systems for which a competitive model can be constructed so that their 
structural stability properties are implied by the same properties of the model. Such a construction 
is possible in the context of the hierarchic stability analysis, and we will show how the decomposition- 
aggregation scheme (refs. 11-14) can be used to form a competitive aggregate model for a large 
class of dynamic systems. The scheme is based on the modern mathematical machinery of the com- 
parison principle and vector Liapunov functions (ref. 15). The decomposition-aggregation stability 
analysis not only provides a competitive aggregate model with structural stability properties, but it 
can also take advantage of the special structural features of complex dynamic systems and yield 
considerable simplification in the relevant stability investigations. 

Let us immediately recall the fact stated in reference 7 that a natural generalization of the 
competitive models considered in reference 3 is represented by a differential equation 

z = u(t,z) (27) 
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where the function u : T X R s R s is u(t, x) G C^°>°) (T X R n ) and u(/, 0) = 0, Vt G 7\ The func- 
tion u(t, x) belongs to the following class of functions: 

K: up, a) < up, b) , Vi - \,2, . . . ,s 

V\{t, a), «,b)\GTXR s B dj = b t , a j <b j (j=\,2,...,s,i*j). 

The class of functions K, used by Kamke (ref. 25) in a formulation of the comparison principle, 
plays an important role in the stability analysis of dynamic systems by the vector Liapunov function. 
Therefore, the strong stability results obtained in the theory of differential inequalities and vector 
Liapunov functions can now be used to study stability properties of the competitive models. 

In the following development, we will show how a competitive model described by equa- 
tion (27) can be constructed for a dynamic system, 

x=f(t,x). (28) 

In (28), x(t) G R n is the state vector of the system and the function/: T X R n -*■ R n is defined, 
bounded, and continuous on the domain T X R n so that the solution x(t\ t 0 , x 0 ) of equation (28) 
exists for all initial conditions (t 0 , x 0 ) G T X R n and t G r 0 . Furthermore, we assume that 
f(t, 0) = 0, Vt G T, and x = 0 is the unique equilibrium state of the system (28). We derive condi- 
tions under which stability of the trivial solution z = 0 of equation (27) implies connective stability 
of the equilibrium x = 0 of the system (28). 

To introduce the connective aspect of stability in the context of system (28), let us represent 
the state space R n as 

R n = R* 1 XR n 2 X ... X R ns (29) 

so that the state vector has the form 

x(t) = [x, T (t) x 2 T (t) . . .x s T (t)] T (30) 

and x ; (r) G R n K The function /(/, x) is further specialized by its components as 

f i (t,x)=f i (t,x i ,e il x 1 ,e i2 x 2 ,. . . ,e is x s ) , / = 1,2, ...,s (31) 

where (r) are elements of the s X s interconnection matrix E(t). Now, we can broaden the scope 
of definition 1 to include the class of systems (28). We say that the equilibrium x = 0 of the 
system (28) is connectively stable if and only if it is stable (in the sense of Liapunov) for all inter- 
connection matrices E(t) (ref. 12). 

To establish connective stability of the system, we will form the model (eq. (27)) where the 
“aggregate” function u{t, x) has the form 

up, z)=up, z u e il z i , e h z 2 , . . . , e (s z s ) , i = 1, 2, . . . ,s (32) 

and Zj- G R 1 . On the basis of (32), we define the model 
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z = u(t, z ) 


(33) 


where the function n(r, z) 6 K corresponds to the fundamental interconnection matrix E. Then, 
stability of the trivial solution z = 0 of the aggregate competitive model (33) implies connective 
stability of the equilibrium x = 0 of the system equation (28). 

More precisely, we prove the following: 

Theorem 3. There exists a function v(t, x) with the properties: 

v(t, x) E C^°’° ) (T X R n ); v(t, x) is locally Lipschitzian in x; v(t, 0) = 0, and 

v(t, x)> 0 on T X R n ; the function v: T X R n -*• R + l , defined as 

s 

v{t,x) = b T v(t,x ) = ^2 b i v i( r ’ x ) ( 34 ) 

i=\ 

for some constant vector b > 0 satisfies the inequalities 

0/(lkll)<Kt,x)<0 / /||x||), V(t,x)ETXR n (35) 

where 0/, 4>[j e H and (pjip) -*■ + 00 as p -*■ + the function 

D + v(t,x) = lim sup f {v[r +/?,x + ft/(/,x)] -v(r r x)[ (36) 

h -K)+ 

defined with respect to equation (28), satisfies a differential inequality 

D + v(fx)^U[t, v(tpc)] , V(t,v)ETXR s (37) 

vv/iere the function u: T X R + s R s is u(t, v) E K and u(t, 0) = 0. 

Then, asymptotic stability in the large of z = 0 in the comparison equa- 
tion ( 33) implies connective asymptotic stability in the large of x = 0 in 
equation ( 28), and 


v(t,x) = [v,(t, x) v 2 (t,x) . . . v s (r,x)] T (38) 

is a vector Liapunov function for the system (28). 

Proof : It is a well-known result (ref. 1 5) that under the conditions of the theorem, asymptotic 
stability of z = 0 of equation (33) implies the same stability property forx = 0 in equation (28). 

To show that stability of x = 0 is also connective, we notice that the assumption v(t , x) > 0 on 
T X R n allows us to establish the inequality 

u(t, z) < u(t, z ) , V(t, z)ET X R s (39) 


which holds for all Eft). Therefore, the differential inequality in equations (37) is also valid for all 
interconnection matrices £(/), and stability of x = 0 is connective. This proves theorem 3. 
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An attractive property of theorem 3 is that it allows a reduction in dimensionality of stability 
problems in much the same way as does the original Liapunov theory. By theorem 3, we can prove 
stability of a system of order n by demonstrating stability of s lower-order systems and one system 
of order s < n which involves the vector Liapunov function. 

Another important feature of theorem 3 we will not take advantage of is that the components 
Vj(t, x) of the vector Liapunov function v(t, x) can satisfy the weaker requirements than the usual 
ones associated with scalar Liapunov functions. This is not to say, however, that vector Liapunov 
functions are automatically more flexible than scalar Liapunov functions as was demonstrated in 
reference 7. 


Let us now show how to construct the competitive aggregate model (27) for a given system (28) 
by applying the decomposition-aggregation method proposed in reference 13. We start with decom- 
posing the system (28) into 5 subsystems described by 

x ( = g f (t, x ( ) + hjit, x) , i = 1 , 2, . . . , s (40) 

where xfit) E R n i is the state vector of the subsystem (40) and represents the z'th “vector” component 
of the state vector x(r) specified by equation (30). 

In equation (40), the functions gj : T X R ni -* R n ‘ describe the “decoupled” subsystems: 

x i =g i (t,x i ), i=l,2, ...,s (41) 

obtained from equation (40) when all interaction functions : T X R n -*■ R ni among the subsys- 
tems are set to zero (hjfj, x) = 0). Each g^t, 0) =0, Vt E T, so that x ; - = 0 is the unique equilibrium 
state of every subsystem (41). 


We assume that with each decoupled subsystem (41) we can associate a function 
Vj : T X R n i -* R\. such that v^t, xf) E c( 0>0 ) ( T X /?”*), Vj(t, x) satisfies a Lipschits condition in 
Xj for a constant > 0, and 

0/, (Ik/ll) < Vj(t, Xj) < 0/2 (liXjll) 

D\<t, x / )(4i)<^0 /3 (llx / ll) ■ } (42) 

Vi = 1, 2, . . . ,s ; Vt ST , Vx t E R n ‘ 

where D+vfit, , ) = lim sup ^ + h, x ^ + hgfit, x^)] - vfit, x z )}, the functions 0,-, , 0 i2 , 

. h-* 0 + 

0 ; - 3 E H and 0^, (p) -> + °° as p-> + °°. 

We assume that the interactions hj(t, x) among the subsystems (41) have the form 

hj(t, x) = hjit, e if Xi , e h x 2 , . . . , e [s x s ) (43) 

and that there exist bounded functions £ : T X R n -*■ /?_j_ such that 
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( 44 ) 


s 

llfift, *)il < ^2 e ( j{ t) £(/, x) <t> j3 {\\Xj\\) , 

7-1 

Vi= 1,2, . . . ,s ) V(t,x)eTXR n 

> 

Let us define am Xs constant matrix A = (dp) by 


^ij + i K i a ij 


(45) 


where are the elements of the fundamental interconnection matrix E, and the numbers ay.- > 0 
are computed as 


oijj = max 




sup x)}- 

TXR n 


(46) 


Now we state: 

Theorem 4. The equilibrium state x — 0 of the system (28) is connectively 
asymptotically stable in the large if the s X s constant matrix A — ( djj ) defined 
by equation (45) satisfies the conditions (6). 

n 1 

Proof : As in reference 13, consider the function v : T X R -*■ R + , 

s 

v(t, X) = ^2 bjVj(t, Xj) (47) 

i= 1 

as a candidate for a Liapunov function for the system (28) where bj > 0, i = 1, 2, . . . , s are com- 
ponents of a positive (yet unspecified) vector b. Then define the function 

D + v i (t, x z )( 40 ) = lim sup ^ {v ; (t + h, x f + h [gfi, x,-) + hj(t, x)] - vft, x,-)J (48) 

h (j 

which we compute with respect to subsystem equation (40) and obtain 

D + v i (t, x,)( 4 o ) < D+vft, x,)( 4 , ) + K^lhft, x)|| , i = 1,2, ... ,s (49) 

since each vft, x,-) is Lipschitzian with a constant Kj > 0 (ref. 1 5). 

By using the interconnection conditions (44) and subsystem stability equations (42), from 
(49) we obtain the vector differential inequality: 

Z) + v( 40 ) < A w(v) , V(t, x)(=TXR n (50) 

which holds for all interconnection matrices E(t). 
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As in equation (13), w : R s -> R s is the comparison vector function: 

w(v) = [ 0 , 3 ^ 1 ) 023 (^ 2 ) • • • ^ 3(^)1 T ■ (51) 

By using the same argument as in the proof of theorem 1, from equations (47) and (50), we 
obtain 

0/(ll*ll)<KL*)<0//(IMI), D + v(t,x) <-0 //7 dk||), V(t,x)£TXR n (52) 

where the functions 

s 

0/01*11) = b m ^2 0/i (il*/ll) , 

1=1 

s 

0//OWI) = % £ 0 /2 (|l*/ll), 

/= 1 

s 

<P/ni\\x\\) = c m ^2 0/ 3 (ll*/ll) 

/=! 


all belong to the class H. 

Since inequalities (52) hold for all interconnection matrices E(t), the equilibrium state x = 0 
of the system (28) is connectively asymptotically stable in the large. The proof of theorem 4 is 
complete. 

By observing that riw(v) G K, theorem 4 follows. directly from inequality (50) and theorem 3. 
In the proof of theorem 4, we used the scalar function v = b J v to show explicitly the multilevel 
nature of the analysis. The components of the vector Liapunov function v are scalar Liapunov 
functions responsible for stability on the subsystem level, and the scalar function v is a scalar 
Liapunov on the aggregate overall system level. 

If we strengthen the constraints imposed on the interactions hj(t, *) among the subsystems, 
we can infer exponential stability of the overall system from exponential stability of the subsystems. 


Let us assume that the estimates (42) have the form 


1 1 1*/ 1 1 <»/(L */)<t? ( 


12 


| !*,•! | D V t {t, X ; )(3 9 ) < -7 TfVjit, Xj) , 


Vi = 1, 2, . . . , j ; Vt^T, VXi&lf' 


(53) 


where t? ;i , Tfy 2 , and ^ are all positive numbers. Inequalities (53) guarantee exponential stability of 
the decoupled subsystems (39). 
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We also assume that the interactions hf_t, x) among the subsystems satisfy the constraints 

s 

11/2/0, X)\\ < 5T *)ll*/ll Vi=\,2,...,s\ V<f.x)ETXR n . (54) 

7-1 

We again construct the aggregate s X s constant matrix A - (flu) defined by 

a ij ~ ^ij n i (55) 

and prove the following: 

Theorem 5. The equilibrium state x = 0 of the system (28) is connec lively 
exponentially stable in the large if the s X s constant matrix A — (a,y) defined 
by equation ( 55) satisfies the inequalities ( 6). 

Proof. We again use the function v(t, x ) of equation (47) as a Liapunov function. By using the 
estimates (53) and (54), from (49) we obtain the following inequalities: 


D + v i (t, X/)( 40 )^~ir i v i (t, Xj) + K t Y 

7-1 


e if x ij\\x\\ 


s 

<-*?t<t. x i) + Y 

7=1 

i = 1,2, ... ,s ', VtET, 


eifHj-npft, Xj) 

VxER n 


that are valid for all interconnection matrices E(t). 


(56) 


From equation (56), we obtain the aggregate competitive model: 

Z) + v (40 ) <Av , V(t, x) E T X R n (57) 

which is a linear differential inequality in v. Proceeding as in the proof of theorem 4, we derive the 
scalar inequality 


D + v(t, jc) (4 0) < -7 ru(t, x) , V(t. x)ET X R n , V£(t) 


(58) 


where it — min 7 r,- . 

i 1 

Integrating inequality (58), we obtain 


u[t, x(t)] <K/o,^o)exp[-7r(/-7 0 )] , Vt E T 0 , V(t 0 , x 0 ) E T X R n , VE(t) . (59) 
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By use of the estimates (53), we can obtain from this inequality a further inequality that 
involves the solution x(t; t 0 ,x 0 ) of the original overall system: 

lk(t;t 0 . Jc o)IKnilJc 0 llexp[-7r(t-f 0 )] , VteT 0 , V(t 0 , x 0 ) E T X R n 

where n = s 1/2 b m b M r }^ , r\ Ml , b m = inin b ( , b M = rnax b t , r? m , = m.in r?,-, , 
and 7r = min 7r ; - . This proves theorem 5. 


, VE{t) (60) 
2 = “f* ^2 - 


CONCLUSIONS 


Two important conclusions result from this work. First, the competitive structures are an 
appropriate framework for constructing reliable complex systems. Secondly, via the diagonal 
dominance of the aggregate competitive model, the decomposition-aggregation method provides a 
good measure of complexity for stable, large-scale systems. Both results arose in the context of 
connective stability by applying the mathematical apparatus of the comparison principle and vector 
Liapunov function. One of the most attractive aspects of the results is that they are obtained in the 
Hicks-Metzler algebraic setting, which provides a rich environment for their application and further 
refinements. 

A number of interesting problems were either not mentioned or were not explored in sufficient 
detail. For example, it is not clear how the decomposition and aggregation should be performed to 
balance the gains in simplification against the errors resulting from the approximation involved in 
the decomposition-aggregation process. Furthermore, it is of interest to investigate various implica- 
tions (ref. 26) of the obtained results in the competitive analysis conducted in the field of economics 
(refs. 2,3,1) and ecosystems (refs. 4, 8, 10). In control systems (ref. 27), there are already some 
specific applications that indicate certain definite advantages of the competitive analysis and the 
decomposition-aggregation method in the multilevel stabilization and optimization of large-scale 
systems. A good deal of work remains to be done to explore the possibilities offered by the competi- 
tive structures in dynamic systems and to obtain important new results. 
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